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Abstract 

For a fading Gaussian multiple access channel with user cooperation, we obtain the 
optimal power allocation policies that maximize the rates achievable by block Markov 
superposition coding. The optimal policies result in a coding scheme that is simpler 
than the one for a general multiple access channel with generalized feedback. This 
simpler coding scheme also leads to the possibility of formulating an otherwise non- 
concave optimization problem as a concave one. Using the channel state information 
at the transmitters to adapt the powers, we demonstrate significant gains over the 
achievable rates for existing cooperative systems. 
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1 Introduction 



Increasing demand for higher rates in wireless communication systems have recently triggered 
major research efforts to characterize the capacities of such systems. The wireless medium 
brings along its unique challenges such as fading and multiuser interference, which make 
the analysis of the communication systems more complicated. On the other hand, the same 
challenging properties of such systems are what give rise to the concepts such as diversity, 
over- heard information, etc., which can be carefully exploited to the advantage of the network 
capacity. 

In the early 1980s, several problems which form a basis for the idea of user cooperation in 
wireless networks were solved. First, the case of a two user multiple access channel (MAC) 
where both users have access to the channel output was considered by Cover and Leung [1], 
and an achievable rate region was obtained for this channel. Willems and van der Meulen 
then demonstrated [2] that the same rate region is achievable if there is a feedback link to 
only one of the tansmitters from the channel output. 

The capacity region of the MAC with partially cooperating encoders was obtained by 
Willems in [3]. In this setting, the encoders are assumed to be connected by finite capacity 
communication links, which allow the cooperation. Willems and van der Meulen also consid- 
ered a limiting case of cooperation where the encoders "crib" from each other, that is, they 
learn each others' codewords before the next transmission [4]. Several scenarios regarding 
which encoder(s) crib, and how much of the codewords the encoders learn, are treated and 
the capacity region for each case is obtained in [4]. The capacity of such channels are an 
upper bound to the rates achievable by cooperative schemes, since in the case of cribbing 
encoders, the sharing of information comes for free, i.e., the transmitters do not allocate any 
resources such as powers, to establish common information. 

An achievable rate region for a MAC with generalized feedback was found in [5]. This 
channel model is worth special attention as far as the wireless channels are concerned, since 
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it models the over-heard information by the transmitters. In particular, for a two user MAC 
with generalized feedback described by [X\ x X 2 , P(y, Hi, 2/2 l^i, ^2), y x x ^2), where user 1 
has access to channel output Y 1 and user two has access to channel output Y 2 , an achievable 
rate region is obtained by using a superposition block Markov encoding scheme, together 
with backward decoding, where the receiver waits to receive all B blocks of codewords before 
decoding. 

Recently, Sendonaris, Erkip and Aazhang have successfully employed the results of these 
rather general problems, particularly that of generalized feedback, to a Gaussian MAC in the 
presence of fading, leading to user cooperation diversity and higher rates [6] . In this setting, 
both the receiver and the transmitters receive noisy versions of the transmitted messages, 
and slightly modifying the basic relay channel case, the transmitters form their codewords 
not only based on their own information, but also on the information they have received 
from each other. It is assumed in [6] that channel state information for each link is known 
to the corresponding receiver on that link, and also phase of the channel state needs to be 
known at the transmitters in order to obtain a coherent combining gain. The achievable 
rate region is shown to improve significantly over the capacity region of the MAC with non- 
cooperating transmitters, especially when the channel between the two users is relatively 
good on average. 

There has also been some recent work on user cooperation systems under various as- 
sumptions on the available channel state information, and the level of cooperation among 
the users. Laneman, Tse and Wornell [7] have characterized the outage probability behavior 
for a system where the users are allowed to cooperate only in half-duplex mode, and where 
no channel state information is available at the transmitters. For the relay channel, which 
is a special one-sided case of user cooperation, Host-Madsen and Zhang [8] have solved for 
power allocation policies that optimize some upper and lower bounds on the ergodic capac- 
ity when perfect channel state information is available at the transmitters and the receiver. 
For a user cooperation system with finite capacity cooperation links, Erkip [9] has proposed 
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a suboptimal solution to the problem of maximizing the sum rate in the presence of full 
channel state information, where it was also noted that the resulting optimization problem 
is non-convex. 

In this paper, we consider a two user fading cooperative Gaussian MAC with complete 
channel state information at the transmitters and the receiver, and average power constraints 
on the transmit powers. Note that, this requires only a small quantity of additional feedback, 
namely the amplitude information on the forward links, over the systems requiring coherent 
combining [6]. In this case, the transmitters can adapt their coding strategies as a function of 
the channel states, by adjusting their transmit powers [10-12]. We characterize the optimal 
power allocation policies which maximize the set of ergodic rates achievable by block Markov 
superposition coding. To this end, we first prove that the seemingly non-concave optimization 
problem of maximizing the achievable rates can be reduced to a concave problem, by noting 
that some of the transmit power levels are essentially zero at every channel state, which also 
reduces the dimensionality of the problem. By this, we also show that the block Markov 
superposition coding strategy proposed in [5] and employed in [6] for a Gaussian channel 
can be simplified considerably by making use of the channel state information. Due to the 
non-differentiable nature of the objective function, we use sub-gradient methods to obtain 
the optimal power distributions that maximize the achievable rates, and we provide the 
corresponding achievable rate regions for various fading distributions. We demonstrate that 
controlling the transmit powers in conjunction with user cooperation provides significant 
gains over the existing rate regions for cooperative systems. 



2 System Model 

We consider a two user fading Gaussian MAC, where both the receiver and the transmitters 
receive noisy versions of the transmitted messages, as illustrated in Figure H The system is 
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modelled by, 



Yo — v hioXi + a/ /i2o-^2 + 



(1) 



Y\ — \/ ^21X2 + Z\ 



(2) 



Y 2 = \/h X2 X x + Z 2 



(3) 



where Xj is the symbol transmitted by node i, Yi is the symbol received at node i, and 
the receiver is denoted by i = 0; Zi is the zero-mean additive white Gaussian noise at 
node i, having variance of, and y/hij are the random fading coefficients, the instantaneous 
realizations of which are assumed to be known by both the transmitters and the receiver. 
We assume that the channel variation is slow enough so that the fading parameters can be 
tracked accurately at the transmitters, yet fast enough to ensure that the long term ergodic 
properties of the channel are observed within the blocks of transmission [13]. 

The transmitters are capable of making decoding decisions based on the signals they 
receive and thus can form their transmitted codewords not only based on their own infor- 
mation, but also based on the information they have received from each other. This channel 
model is a special case of the MAC with generalized feedback [5] . The achievable rate region 
is obtained by using a superposition block Markov encoding scheme, together with backward 
decoding, where the receiver waits to receive all B blocks of codewords before decoding. For 
the Gaussian case, the superposition block Markov encoding is realized as follows [6]: the 
transmitters allocate some of their powers to establish some common information in every 
block, and in the next block, they coherently combine part of their transmitted codewords. 
In the presence of channel state information, by suitably modifying the coding scheme given 
by [6] to accommodate for channel adaptive coding strategies, the encoding is performed by 




(4) 
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for i, j e {1, 2}, i ^ j, where Ajo carries the fresh information intended for the receiver, 
Xij carries the information intended for transmitter j for cooperation in the next block, and 
Ui is the common information sent by both transmitters for the resolution of the remaining 
uncertainty from the previous block, all chosen from unit-power Gaussian distributions. All 
the transmit power is therefore captured by the power levels associated with each component, 
i.e., Pio(h), Pij(h) and ^(h), which are required to satisfy the average power constraints, 



E\p i0 {h) + Pij {h) + PUl (h)) = E\pi{h)\ < Vi 



1,2. 



(5) 



Following the results in [6], it can be shown that the achievable rate region is given by 
the convex hull of all rate pairs satisfying 
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(6) 
(7) 



where the convex hull is taken over all power allocation policies that satisfy 

For a given power allocation, the rate region in ©-© is either a pentagon or a triangle, 
since, unlike the traditional MAC, the sum rate constraint in (jHJ) may dominate the individual 
rate constraints completely. The achievable rate region may alternatively be represented as 
the convex hull of the union of all such regions. Our goal is to find the power allocation 
policies that maximize the rate tuples on the rate region boundary. 
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3 Structure of the Sum Rate and the Optimal Policies 

We first consider the problem of optimizing the sum rate of the system, as it will also shed 
some light onto the optimization of an arbitrary point on the rate region boundary. The sum 
rate (JBJ) is not a concave function of the vector of variables p(h) = [pio(h) pi2(h) pu x (h) P2o(h) 
P2i(h) p[/ 2 (h)], due to the existence of variables in the denominators. In what follows, we 
show that for the sum rate to be maximized, for every given h, at least two of the four 
components of [pio(h) P12Q1) P2o(h) ]?2i(h)] should be equal to zero, which reduces the 
dimensionality of the problem and yields a concave optimization problem. 

Proposition 1 Let the effective channel gains normalized by the noise powers be defined as 
Sij = hij/o-j. Then, for the power control policy p*(h) that maximizes (QJ) ; we need 

1. p* (h) = p 20 (h) = 0, if s l2 > sio and s 21 > s 20 

2. p* w (h) = p 21 (h) = 0, if s 12 > sio and s 21 < s 20 

3. p* 12 (h) = p 20 (h) = 0, if s 12 < sio and s 21 > s 20 
I p* 12 (h) = p* 21 (h) = 

OR 

Ko(h)=P2i(h)=0 
OR 

p* 12 (h) = p* 20 (h) = r 

Proof: To simplify the notation, let us drop the dependence of the powers on the channel 
states, whenever such dependence is obvious from the context. Let pi = p^ + pij + py i be 
the total power allocated to a given channel state. Let us define 



A = 1 + siopi + s 20 p 2 + 2 v /si S2oPc/ 1 Pc/ 2 (9) 

B = l + ^io + ^ 2 o (1Q) 
(1 + s 12 pio)(l + S21P20) 

C = (1 + S12O10+P12)) (1 + s 21 {p 20 + p 21 )) (11) 



> ifsu < sio and s 21 < s 20 
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Then, an equivalent representation of the sum rate in (JHJ) is 



R snm = min {E [\og{A)\ , E [log(SC)]} (12) 

Now, let us arbitrarily fix the total power level, pi, as well as the power level used 
for cooperation signals, pjj. < Pi, allocated to a given state for each user. For each such 
allocation, the quantities A and C appearing in the sum rate expression are fixed, i.e., 
allocating the remaining available power pi — pu i among p i0 and Pij will not alter these 
quantities. Note that, such allocation also does not alter the total power consumption at the 
given state, so we may limit our attention to the maximization, 

max B(p w ,p 2 o) 

{P10,P20} 

s.t. pio +p 12 = pi -p Ux 

P20 + P21 = P2 - Pu 2 (13) 



The partial derivatives of B with respect to p\o and P20 are 

9B _ sip - si 2 (l + S20P20) , 14 x 

dp 10 (1 + si 2 pio) 2 (l + S21P20) 

9B s 20 - s 2 i(l + gioPio) , . 

dp 20 (1 + s 2 ip 2 o) 2 (l + siaPio) 

Therefore, we make the following conclusions: 

1. S12 > sio, s 2 i > s 20 . Then, g < and g < 0, i.e., S(p 10 ,p 2 o) is monotonically 
decreasing in both p w and p 2 o, therefore the sum rate is maximized at p w = p 20 = 0. 

2. S12 > sio, S21 < S2o- Then, < 0, and the function is maximized at pio = for 
any P20- But this gives | Plo =o > 0, meaning p 2 o should take its maximum possible 
value, i.e., p 2 i = 0. 



8 



3. S12 < sio, S21 > S20 Follows the same lines of case 2 with roles of users 1 and 2 
reversed. 

4. S12 < sio, S21 < S20- In this case, the partial derivatives of .£> can be both made 
equal to zero within the constraint set, yielding a critical point. However, using higher 
order tests, it is possible to show that this solution corresponds to a saddle point, and 
B is again maximized at one of the boundaries, pio = 0, P20 = 0, pio = Pi — Pui, 
P20 = P2 — Pu 2 - Inspection of the gradient on these boundary points yields one of the 
three corner points {{p\ —pu-L, 0), (0,p2 — pu 2 )i (Pi ~PUi,P2 —Pu 2 )} as candidates, each 
of which corresponds to one of the solutions in case 4. 

Although two of the components of the power vector are guaranteed to be equal to 
zero, which ones will be zero depends on the pi and pu i that we fixed, therefore we 
are not able to completely specify the solution, independent of Pi and pu v in this case. 
On the other hand, the settings of interest to us are those where the channels between 
the cooperating users are on average much better than their direct links, since it is 
in these settings when cooperative diversity yields high capacity gains [6]. In such 
scenarios, the probability of both users' direct link gains exceeding their corresponding 
cooperation link gains (case 4) is a very low probability event. Therefore, which of the 
three possible operating points is chosen is not of practical importance, and we can 
safely fix the power allocation policy to one of them to carry on with our optimization 
problem for the other variables. Although admittedly this argument is likely to cause 
some suboptimality in our scheme, as will be seen in the numerical examples, we still 
obtain a significant gain in the achievable rates. 

□ 

The significance of this result is two-fold. Firstly, given a channel state, it greatly sim- 
plifies the well known block Markov coding, in a very intuitive way: if the direct links of 
both users are inferior to their cooperation links, the users do not transmit direct messages 
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to the receiver as a part of their codewords, and they use each other as relays. If one of 
the users' direct channel is better than its cooperation channel, and the other user is in the 
opposite situation, then the user with the strong direct channel chooses to transmit directly 
to the receiver, while the weaker direct channel user still chooses to relay its information 
over its partner. Second important implication of this result is that it now makes the prob- 
lem of solving for the optimal power allocation policy more tractable, since it simplifies the 
constraint set on the variables, and more importantly, this makes the sum rate a concave 
function over the reduced set of constraints and variables. 

Corollary 1 The sum, rate R sum given by Ufy . is a conca/ve junction o/p(h) ; over the 
reduced constraint set described by Proposition 1. 

Proof: The proof of this result follows from directly substituting the zero power components 
into the sum rate expression in (jl2j) . Note that in each of the four cases, the second function 
in the minimization, i.e., log(-BC) takes either the form log(l + a) + log(l + 6), or log (1 + a +6), 
both of which are clearly jointly concave in a and b. Also, log (A) is clearly a concave function 
of p(h) since it is a composition of a concave function with the concave and increasing 
logarithm. The desired result is obtained by noting that the minimum of two concave 
functions is concave. □ 

Thus far we have discussed the structure of the sum rate, as well as some properties of 
the optimal power allocation that maximizes that rate. We now turn back to the problem of 
maximizing other rate points on the rate region boundary. To this end, we point out another 
remarkable property of the solution in Proposition 1. Consider maximizing the bound on 
R\ in (jnj). For fixed pu 1 and pi, it is easy to verify that all of the available power should be 
allocated to the channel with the higher gain, i.e., if > Sio, then we need pxo = an d 
Pn = Pi —Pu-i- The same result also applies to R 2 . But this shows that, the policies described 
in Proposition 1 completely agree with optimal policies for maximizing the individual rate 
constraints in cases 1-3, and they also agree if we choose the operating point in case 4 to 
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be P12 = 0, P21 = 0. Therefore, the allocation in Proposition 1 enlarges the entire rate 
region in all directions (except for the subtlety in case 4 for the sum rate). This has the 
benefit that the weighted sum of rates, say _R M = HiR\ + /i 2 -R 2 also has the same concavity 
properties of the sum rate, since for /ij > fij, the weighted sum of rates can be written as 
Rfj. = /^j-Rsum + (yUi — Hj)Ri, where both -R S um and Ri are concave. Optimum power control 
policies that achieve the points on the boundary of the achievable rate region can then be 
obtained by maximizing the weighted sum of rates, which is the goal of the next section. 



4 Rate Maximization via Subgradient Methods 

In this section we focus on maximizing the weighted sum of rates. To illustrate both the 
results of the preceding section and the problem statement for this section more precisely, 
let us consider, without loss of generality, the case when fix > /i 2 , and write down the 
optimization problem explicitly: 

max (/i x - fi 2 ) \ Ei 2 [log(l +pi 2 (h)si 2 )] 

P(h) [ 

+ E 3A pog(l+pio(h)sio)]J+^niin|s[log(A)], 

+ E 1 [log(l +p 12 (h)s 12 ) + log(l +P2i(h)s 2 i)] 
+ E 2 [log(l + Pi 2 (h)s 12 ) + log(l + p 20 (h)s 20 )] 
+ E 3 [log(l +pio(h)si ) + log(l +j?2i(h)s 2 i)] 
+ Ei [log(l +pio(h)si +P2o(h)s 20 )] I 

s.t. £3,4 \pw(h)) + E lj2 [pi 2 (h)] + E [puAh)} < p l 

E 2A [p2o(h)] + E li3 [p 2 i(h)] + E [ PU2 (h)} < p 2 (16) 

where, Es denotes the expectation over the event that case S C {1, 2, 3, 4} from Proposition 
1 occurs, and A is as given by Note that the objective function is concave, and the 
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constraint set is convex, therefore we conclude that any local optimum for the constrained 
optimization problem is a global optimum. Although is different iable almost everywhere 
since it is concave, its optimal value is attained along the discontinuity of its gradient, namely 
when the two arguments of the min operation are equal. Hence, we solve the optimization 
problem using the method of subgradients from non-differentiable optimization theory [14, 
15]. 

The subgradient methods are very similar to gradient ascent methods in that whenever 
the function is different iable (in our case almost everywhere), the subgradient is equal to the 
gradient. However, their major difference from gradient ascent methods is that they are not 
necessarily monotonically non- decreasing. Let us denote the objective function in (JTfij) by 
A subgradient for the concave function -R M (p) is any vector g that satisfies [15], 

RAp') < R M + (p' - p) t s (17) 

and the subgradient method for the constrained maximization uses the update [15] 

p(k + l) = [p(k) + a k g k }i (18) 

where denotes the Euclidean projection onto the constraint set, and a k is the step size 
at iteration k. There are various ways to choose a k to guarantee convergence to the global 
optimum; for our particular problem, we choose the diminishing stepsize, normalized by the 
norm of the subgradient to ensure the convergence [14] 

a 1 . 

While the subgradient method is a powerful tool to obtain the global optimum of (fTHjl. it 
does not shed light onto the characteristics of the optimal power allocation policies. Although 
it is in general difficult, if not impossible, to state the optimal policy in a closed form, we 
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provide a simplified case as an example, and provide conditions the optimal policies should 
satisfy, in the next subsection. 

4.1 Properties of the Optimal Power Allocation: An Example 

Throughout the analysis in this subsection, we assume for simplicity of the derivations that 
the fading distributions are such that all realizations of the fading values satisfy si 2 > s±o 
and S21 > S20, i-e., the system always operates in the more interesting case 1. This particular 
case is of practical interest since the cooperating transmitters are likely to be closely located 
with less number of scatterers and obstructions when compared to their paths to the receiver, 
and thus have better channel conditions among each other. We also consider maximizing the 
sum capacity, i.e., /ii = /x 2 - Under these assumptions, since Pio(h) = £>2o(h) = always, the 
problem of maximizing the sum rate, given in (JHJ) and also obtained from (j!6)) by inserting 
/ii = /i 2 , reduces to 

max min{£' [log(D)] , E [log (1 + si 2 pi 2 (h)) + log (1 + s 2 iP2i(h))]} 

P(h) 

s.t. E\p 12 (h)+ PUl (h)} <p, 
E\p 21 {h)+p U2 {h))<p 2 
p 12 (h), m (h),p 2 i(h),R/ 2 (h) > 0, Vh 

where 

D = 1 + sio (pi2(h) +PuM) + s 20 (p 2 i(h) +puM) + 2^s 10 s 2oPUl (h)p U2 (h) (21) 

Note that, the maximum of the objective function is attained only when the arguments 
of the min operation are equal, or else it would be possible to increase the smaller argument 
while decreasing the larger, just by transferring the transmit power between the cooperation 
signals Pu l O a ) an d Py(h) at some h, while keeping their sum constant, and therefore not 
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(20) 



violating the average power constraints. Hence, to overcome the difficulties arising from the 
non-differentiability of the cost function, the min in the cost function can be replaced by 
an equality constraint in the constraint set, and the problem can be restated with a new 
differentiable objective function as 



max E [log (1 + s 12 p 12 (h)) + log (1 + s 21 p 21 (h))] (22) 

P(h) 

s.t. E [log(D)] = E [log (1 + s 12 p 12 (h)) + log (1 + s 21 p 21 (h))] (23) 

E\p l2 0a)+p Ul {h)]<p 1 (24) 

E[p 21 (h)+p U2 (h)]<p 2 (25) 

p 12 (h), PUl (h),p 21 (h),p U2 (h) > 0, Vh (26) 



Associating the Lagrange multiplier 7 to the equality constraint (|23j) . Ai,A 2 > to the 
power constraints (|24 |) . (|25j ) . and £i(h) > 0, i — 1, • • • , 4 to the non-negativity constraints 
(|2f)j) . and noting that the power constraints need to be satisfied by equality, we obtain the 
Karush-Kuhn- Tucker (KKT) necessary conditions for optimality. Note that for optimality, 
for any given h, the components pu x (h) and pu 2 (h) should be either both positive or both 
zero. Therefore, it is sufficient to analyze these two cases. 

First let pu x (h) , pu 2 (h) > 0. Then, the KKT conditions reduce to 



(l + 7) :r- SU ^ <Ar+7f (27) 
l + s 12 p 12 {h) D 

(1 +7) T -^- 7 rv <A 2 + 7^ (28) 

Vgiog2oPi/ 2 (h) + sioy^(h) = _^ 
7 ^y^ (h) 
V / SioS2oPc/ 1 ( h ) + smy/ Pu a (h) = _ x 

From (|29p and (|30|) . it follows that there is a linear relationship between the optimal ^^(h) 
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and pu 2 (h), i.e., 



(31) 



Next, we include the case when pu 1 (h) = pu 2 (h) = in our analysis. The partial 
derivatives g p ^ h ) , i — 1,2 are not well defined when pc/^h) are both zero. However, the 
above derivation shows that the optimal solution is guaranteed to lie on the line given by 
(PU, since (p[/i(h.),pj/ 2 (h)) = (0,0) also lies on the same line. Therefore, it is sufficient to 
evaluate the optimality conditions along the line (J3*T|) for pu^h) = pt/ 2 (h) = 0, as well as for 
p [ / 1 (h),p C / 2 (h) > 0. Substituting flZJ into (|2Z), into using (jSTJ), and extending flZJ) 
and (J3Uj) to include the case p{/ 1 (h),p[/ 2 (h) = (0,0), we obtain the overall set of necessary 
conditions for the optimality of the power allocation policy as, 



Sl2 



< 



1 



A1S20 



1 + Sl2Pl2(h) (1 + 7) A 2 Sio + A1S20 
S21 



< 



1 A|sio 



1 + S 2 lP2l(h) (1 + 7) A 2 Sio + A1S2O 

\ 7 ( 1 Al 

_ < d [ Sw + Y 2 S2 ° 

-A 2 < - ( .20 + ^.10 



w.e. ifpia(h) > (32) 



w.e. if p 2 i(h) > (33) 
w.e. ifpc/^h) > (34) 



w.e. ifpc/ 2 (h) > (35) 



where w.e. stands for "with equality." 

These conditions shed some light onto the structure of the optimal power allocation 
policy: when the channel gains sjo, S20 from the transmitters to the receiver are fixed, the 
optimal power levels associated with the components intended for the other transmitter, 
Pi2(h) and p 2 i(h), are independent of each other and are obtained by single user waterfilling 
over the individual link gains S12 and S21, respectively. The direct link gains sio and S20 
determine the water level. This is easily observed by denoting the right hand sides of f!32l) 
and (JBlIj) by the constants V\ and v 2 (since s w and s 2 o are fixed), and solving for pi 2 (h) and 
p 2 i(h), keeping in mind the equalities are satisfied only when the power levels are positive, 
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i.e., 



Pn(h) =(--—) 
P2i(h) =(--—) 

\V2 S21J 



+ 



(36) 



+ 



(37) 



where (-) + denotes max(0, •). 

Once pu(h) and £>2i(h) are calculated, they can be substituted into and (|33j) with D 
defined as in (|21|) to obtain pjj. . On the other hand, the computation of the optimal power 
levels require solving for the Lagrange multipliers Aj and 7 that satisfy the power constraints, 
which would require a multi-dimensional search. Moreover, due to the constraint (J28)) being 
potentially non-convex for some power values, the solution to the KKT conditions may not be 
sufficient to guarantee optimality. In the next section, we obtain the optimal power allocation 
policy numerically using the subgradient method, and then demonstrate by simulations that 
the optimal power allocation policy possesses the properties described in this subsection. 

5 Simulation Results 

In this section we provide some numerical examples to illustrate the performance of the 
proposed joint power allocation and cooperation scheme. 

Figure El illustrates the achievable rate region we obtain for a system with p { = of = 1, 
subject to uniform fading, where the links from the transmitters to the receiver are sym- 
metric and take values from the set {0.025, 0.050, •• • ,0.25}, each with probability 1/10, 
while the link among the transmitters is also symmetric and uniform, and takes the values 
{0.26, 0.27, • • ■ , 0.35}. Notice that here, we have intentionally chosen the fading coefficients 
such that the cooperation link is always better than the direct links, therefore, the system 
operates only in case 1 of Proposition 1. Consequently, in this particular case, our power 
allocation scheme is actually the optimal power allocation policy for the block Markov su- 
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perposition encoding scheme. 

The region for joint power control and cooperation is generated using the subgradient 
method with parameters a = 50 and 6 = 5. We carried out the optimization for various 
values of the priorities /i, of the users, each of which give a point on the rate region boundary, 
and then we performed a convex hull operation over these points. We observe that power 
control [10-12] by itself improves on the rate region of the cooperative system with no power 
control, for rate pairs close to the sum rate, by utilizing the direct link more efficiently 
Joint user cooperation and power control scheme significantly improves on all other schemes, 
as it takes advantage of both cooperation diversity and time diversity in the system. In 
fact, we can view this joint diversity utilization as adaptively performing coding, medium 
accessing and routing, thereby yielding a cross-layer approach for the design of the wireless 
communication system. 

Figure El also corresponds to a system with unit SNR, but this time subject to Rayleigh 
fading, i.e., the power gains to the receivers are exponential random variables, with E[hio] = 
E[h 2 o] = 0.3, E[hi 2 ] = E[li2i] = 0.6. In this setting, all four cases in Proposition 1 are 
realized, and there is potentially some loss over the optimally achievable rates. However, we 
obtain a very similar set of rate regions to the uniform case, indicating that in fact the loss, if 
any, is very small thanks to the very low probability of both of the direct links outperforming 
the cooperation link. 

It is interesting to note in both Figures 121 and El cooperation with power control improves 
relatively less over power control only near the sum capacity. This can be attributed to the 
fact that, for the traditional MAC, the sum rate is achieved by time division among the 
users, which does not allow for coherent combining gain [11]. Therefore, it is not surprising 
to see that in order to attain cooperative diversity gain, users may have to sacrifice some of 
the gain they obtain from exploiting the time diversity. 

In Figure El we illustrate the convergence of the subgradient method. The objective 
function is with /!% = 2 and fi 2 = 1, and the step size parameters are varied. We observe 
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that, by choosing larger step sizes, the non-monotonic behavior of the subgradient algorithm 
becomes more apparent, however the convergence is significantly faster than the smaller 
step sizes, as the algorithm is more likely to get near the optimal value of the function in 
the initial iterations. Note that, in our simulations we terminated the algorithm after 1000 
iterations, and the three curves would eventually converge after sufficiently large number of 
iterations. 

Figures 5(a)| and |5(b)| demonstrate the waterfilling nature of the optimal power levels 
Pi 2 (h), and p2i(h), as obtained in Section 4.1. These power levels correspond to the sum 
rate point on Figure El and were obtained by using the subgradient algorithm, after 1000 
iterations. The direct link gains are fixed to the values s w = 0.2 and s 2 o = 0.15. As dictated 
by (JH2*|) and the "water level" for p 12 is higher, and the threshold channel state level 
to start transmitting is lower. We see that after 1000 iterations, the subgradient approach 
nearly satisfies the optimality conditions. 



6 Conclusions 

We have addressed the problem of optimal power allocation for a fading cooperative MAC, 
where the transmitters and the receiver have channel state information, and are therefore 
able to adapt their coding and decoding strategies by allocating their resources. We have 
characterized the power control policies that maximize the rates achievable by block Markov 
superposition coding, and proved that, in the presence of channel state information, the 
coding strategy is significantly simplified: given any channel state, for each of the users, 
among the three signal components, i.e., those that are intended for the receiver, for the other 
transmitter, and for cooperation, at least one of the first two should be allocated zero power 
at that channel state. This result also enabled us to formulate the otherwise non-concave 
problem of maximizing the achievable rates as a concave maximization problem. The power 
control policies, which are jointly optimal with block Markov coding, were then obtained 
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using subgradient method for non-differentiable optimization. The resulting achievable rate 
regions for joint power control and cooperation improve significantly on cooperative systems 
without power control, since our joint approach makes use of both cooperative diversity and 
time diversity. 
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Figure 1: Two user fading cooperative MAC. 



21 



0.35 



0.3 



0.25 



0.2 



0.15 



0.1 



0.05 



Cooperation and Power Control in Uniform Fading, h10=[0.025:0.025:0.25] h12=[0.26:0.01 :0.35] 
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Figure 2: Rates achievable by joint power control and user cooperation for uniform fading. 
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Cooperation and Power Control in Rayleigh Fading, E[h,.]=0.3, E[h 12 ]=0.6 
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Figure 3: Rates achievable by joint power control and user cooperation for Rayleigh fading. 
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Convergence of the Subgradient Method for Different Step Size Parameters 
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^ure 4: Convergence of using subgradient method for different step size parameters. 
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Figure 5: Power distributions obtained using the projected subgradient method. 
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